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Correlation Effects in Carbon Nanotubes 
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We consider the effects of Coulomb interactions on single-wall carbon nanotubes using an on-site 
Hubbard interaction, u. For the (N,N) armchair tubes the low energy theory is shown to be identical 
to a 2-chain Hubbard model at half-filling, with an effective interaction tijv = u/N. Umklapp 
scattering leads to gaps in the spectrum of charge and spin excitations which are exponentially small 
for large N. Above the gaps the intrinsic nanotube resistivity due to these scattering processes is 
linear in temperature, as observed experimentally. The presence of "d-wave" superconductivity in 
the 2-chain Hubbard model away from half-filling suggests that doped armchair nanotubes might 
exhibit superconductity with a purely electronic mechanism. 



PACS numbers: 71.10Hf, 71.10Pm, 78.66Tr 

Carbon nanotubes constitute a novel class of quasi one- 
dimensional (Id) materials which .offer the potential for 
both new physics and technology.^ Although built only 
with carbon atoms, they can be grown in a tremendous 
variety of shapes and sizes. The simplest single- wall tube 
consists of a single graphite sheet which is curved into a 
long cylinder, with a diameter which can be smaller than 
lnm. Several groupsa have succeeded in measuring the 
resistance of a single multi-wall nanotube, composed of 
several concentric cylinders. Crystalline "ropes" consist- 
ing of a triangular packing of (nominally) identical single- 
wall tubes are also very promising, exhibiting signatures 
of metallic transport .0 

Generally, single wall tubes can be characterized by 
two integers, (N,M), which specify the super-lattice 
translation vector which WjLaps around the waist of the 
cylinder. Current theorieserO consist of band structure 
calculations and predict a rich variety of behavior, rang- 
ing from metallic "armchair" tubes with (N, N) to insu- 
lating "zig-zag" tubes with (JV, 0). For very small nan- 
otubes, however, electron correlation effects should be- 
come important, as in other Id systems.Q In this paper 
we study these effects using a tight-binding description 
(which correctly reproduces the band-structure calcula- 
tions) supplemented by an on-site Hubbard interaction 
u. For the (N,N) armchair tubes we show that the ef- 
fective description at low energies is identical to a 2-chain 
Hubbard model at half-filling with an effective interaction 
strength, un — u/N. Since this effective interaction is 
weak for N ~ 10, its effects can be treated perturbatively. 
Particularly important are electronic Umklapp scattering 
processes, present at half-filling. These are predicted to 
open a small charge and spin-gap, changing the behav- 
ior from metallic to insulating at low enough tempera- 
tures. Similar conclusions have been reached indepen-, 
dently in very recent work by Krotov, Lee and Louie.Q 
At temperatures above the charge gap, a simple weak- 
coupling analysis of these interactions gives a resistivity 
which varies linearly with temperature, which may ex- 
plain the observed behavioitl in single carbon "ropes" . 
Furthermore, doping an armchair nanotube is equivalent 
to moving away from 1/2-filling in the 2-chain Hubbard 



model. This problem has been extensively studiedj3 and 
exhibits superconducting behavior at low temperatures, 
with a "d-wave" symmetry This suggests a possible elec- 
tronic mechanism for superconductivity in doped nan- 
otubes. 
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FIG. 1. Illustration of the graphite lattice, with labeling 
and periodic boundary conditions for an armchair tube. 

Following various authors ,0 we first consider a single 
sheet of graphite, composed of carbon atoms arranged on 
the sites of a honeycomb lattice. The underlying Bravais 
lattice is triangular, with two sites per unit cell. The two 
primitive Bravais lattice vectors are a± = (a/2)(±l, y/3), 
where a = ^/3d, with d the near-neighbor carbon sepa- 
ration. Of the four outer shell electrons of each atom, 
three form the sp2 bonds of the lattice, while the fourth 
can tunnel between neighboring p z orbitals. A simple de- 
scription, which correctly accounts for the semi-metallic 
behavior of graphite, consists of a tight binding model 
with one p z orbital per carbon, and a tunneling matrix 
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element t between neighboring atoms. The Bloch states 
for this tight-binding model form two bands, with ener- 
gies E±(k) = ±[|£(fc)| 2 ] 1/2 where 



£(jfe) = 2tcos{k x a/2)e lk y a ' 2 ^ + txe^^, 



(1) 



and k is the crystal momentum. Here we have allowed for 
a different hopping strength, t±, in the y-direction (see 
Fig. 1). With one electron per carbon atom, the Fermi 
energy is at E = 0, with the lower band full and the upper 
empty. The striking feature of this band structure is that 
there are two isolated points in the first Brillouin zone, 
denoted K±, where the bands touch E — 0, and there 
are gapless excitations. In the vicinity of these "Dirac" 
points, for q = k — K± small, the dispersion is relativistic, 
with E(q) = v\q\ and v = (V3a/2)t (for t± = t). When 
t = t±, the gapless points occur at K± = (±47r/3a,0), 
but are shifted along the k x axis for t ^ t± (see Fig. 2). 
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FIG. 2. Dirac points in the Brillouin zone. Dark circles 
and crosses indicate the locations of the gapless points for 
t± — t, while gray symbols schematically indicate the shifted 
positions for t± < t. Dashed lines cut the zone at a discrete 
set of allowed transverse momenta in (a) the armchair tube 
and (b) the zig-zag tube (here with N = 3). 

Single wall nantotubes consist of rolling the honey- 
comb sheet of carbon atoms iato a cylinder. Each tube is 
characterized by two integersQ (N, M), which specify the 
super-lattice translation vector T/n,m) — Na + + Ma*-, 
which wraps around the waist of the cylinder. The crys- 
tal momentum transverse to the axis of the cylinder is 
then quantized. Band structure predictaa metallic be- 
havior whenever the gapless points in the Brillouin zone 
lie on the allowed transverse quantized wavevectors. For 
the armchair nanotubes with (N, N) this is illustrated 
in Fig. 2a, where the allowed values of k y are shown 
as dashed lines for N — 4. Since gapless modes are 
present at k y = 0, band structure predicts metallic be- 
havior for armchair tubes, independent of N. Due to 
curvature effects,B the hopping matrix elements along (t) 
and around (t±) the nanotube will differ slightly, by an 
amount of order 1/N 2 . This shifts the Dirac points at K± 
along k x , but leaves the armchair tube gapless (metal- 
lic). For the (N,—N) zig-zag tube (equivalent to the 
(N, 0) tube) with non- integer N/3, the gapless points do 
not coincide with quantized transverse momenta, so that 
insulating behavior is predicted with a gap varying as 



1/N. For integer N/3 the Dirac points for t = t± arc 
at quantized transverse momenta, but are shifted away 
slightly, of order 1/N 2 , due to curvature effects^ (i =/= t±). 
Thus band structure predicts semi-metallic behavior for 
integer N/3 zigzag tubes (Fig. 2b). 
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FIG. 3. One-dimensional spectrum near Dirac points. 

For the armchair tubes the low energy modes occur 
near the two gapless points, at k y = 0. The Id disper- 
sion away from these two points is shown in Fig. 3. In 
addition there are gapped modes at k y ^ 0, with an en- 
ergy of order An = v(2ir / y/3N) 2 . Below Ajv, the mode 
structure is equivalent to a Id 2-band model, independent 
of the nanotube size N. 

Coulomb interactions can be incorporated into the 
nanotube tight binding model, and will introduce inter- 
actions into the effective Id model. To be concrete we 
focus on an on-site Hubbard interaction, u, so that the 
full Hamiltonian becomes, 



H 



E 

<rr'> 



irr'4 (r)ca (O + u ^ n T (r)n L (r) , (2) 



where the first sum is over spin states (a =|, I) and near 
neighbor sites of the honeycomb lattice, and n a (r) = 
cJ e (r)c Q ,(r). We now show that for (N, N) armchair tubes 
the effective interacting Id model is identical to a 2-chain 
Hubbard model with an interaction strength, un = u/N . 
To do so, we choose a particular basis of states spanning 
the space of low energy states with k y = 0: 



(t>n\{x,y) 



N 1/2 Sx,na S yM ao/Va n even 
N- 1/2 S x , nao 5 y (6e+1)ao/V 3 n odd 



(3) 



{N 1 / 2 6 r rwnfi ran n\ / /n ti even 
N i,2 8 x , nao S y {6e+3)ao/V 3 nodd 

where the second Kronecker delta function must be sat- 
isfied for some integer £, and Oo = a/2. As indicated in 
Fig. 1, <f> n i and 4> n 2 are simply the two normalized basis 
states with uniform support at x — nao on even or odd 
chains, respectively. In the low-energy theory, we may 
restrict the expansion of the field operators to this basis: 



(5) 



Inserting this into the Hubbard Hamiltonian and sum- 
ming over y for fixed x gives 



H = 



| ^ ( c liia c n+l,i 



nil ™T ™i n d' 



*-L ( c lla C n2a +^- C -)} 
(6) 
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which is precisely the Hamiltonian of the two-chain Hub- 
bard model, but with an effective weak interaction um — 
u/N . The factor of 1/N arises because the electrons are 
delocalized around the circumference of the nanotube, 
and hence occupy the same site with a probability re- 
duced by 1/N. 

A considerable amount is known about the two-chain 
Hubbard model,E3 particularly in the weak-coupling limit, 
where controlled renormalization group calculations may 
be used. These methods proceed by diagonalizing the 
kinetic energy, linearizing the Id spectrum near the re- 
sulting Fermi points, and expanding the Hubbard inter- 
actions in the basis of states of the resulting two bands. 
In the undoped case at half-filling, interactions drive an 
instability to a Mott-insulating spin liquid with a gap 
in both the charge and spin sectors. In the weak cou- 
pling limit, un <C t, both gaps are exponentially small: 
A ~ texp(—ct/uN)- At temperatures below the charge 
gap, A c , activated behavior is expected in the resistivity, 
p ~ exp(A c //cB? 1 ), as illustrated schematically in Fig. 5. 
With increasing un the charge gap evolves continuously 
into the strong coupling Mott gap associated with the 
energy cost of doubly occupying a site. The spin gap 
at strong coupling is more subtle, but indicates a quan- 
tum disordered or short-range resonating valence bond 
(RVB) ground state.El A spin gap is also present in a 
two-leg Heisenberg ladder, in contrast to a single chain 
which has gapless spin excitations.EH For the armchair 
tubes the spin gap at strong coupling can be understood 
in terms of a Heisenberg spin model on the honeycomb 
network, which is topologically equivalent to the "brick 
wall" lattice shown in Fig. 4a. In the anisotropic limit 
J± ^> J, local spin-singlets form along the vertical rungs, 
and there is a spin-gap to triplet excitations. 
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FIG. 4. Effective spin models for (a) the armchair tube and 
(b) the zig-zag tube. 

Upon doping for t± < 2t, the two-chain Hubbard 
model is known to undergo a phase transition into a state 
which retains the spin gap but develops power-law singlet 
superconducting (SS) and charge-density-wave (CDW) 
correlations. □ Furthermore, the pair wavefunction asso- 
ciated with the SS correlator has approximate d x 2_ y 2 
symmetry (i.e. a sign change from quadrant to quad- 
rant ini-tie k x , k y plane). Both theoretical and numerical 
studiesuO suggest that the SS correlations are enhanced 
over the CDW ones if the Fermi level is pushed into the 
proximity of a band edge. In weak coupling, the enhance- 
ment is mediated by scattering into the nearly empty /full 
band, for which the Id van Hove singularity provides an 



enormous density of states. This suggests that for arm- 
chair tubes, superconducting effects might be maximized 
by tuning the doping so that the Fermi energy coincides 
with the lowest- lying (k y = ±2n/ Nay/3) bands near the 
two Dirac points (see Fig. 3). 

For (JV, — N) zig-zag tubes with integer 7V/3 the gap- 
less Dirac points coincide with the discrete quantized 
momenta for t = t±, see Fig. 2b. However, due to 
curvature effects t / t±, so that the gapless points 
are slightly shifted leading to a small gapjj of order 
\t — ijj ~ 1/N 2 . Since this is smaller than the effective 
interaction strength, which varies as 1/N, it is probably 
legitimate to ignore this small shift when interactions are 
included. Even with this simplification, it is not possible 
to map the zig-zag tube directly into a two-chain Hub- 
bard model. Nevertheless, proceeding by focusing on the 
two gapless modes and expressing the Hubbard interac- 
tion in terms of these, one obtains an effective interacting 
Id two-band model for the zig-zag tube. Just as for the 
two-chain Hubbard model, this model has Umklapp scat- 
tering processes, but their particular strengths are differ- 
ent. It is natural to expect that these Umklapp processes 
will gap out both the charge and spin excitations, just as 
for the Hubbard model, although a definitive statement 
requires a detailed calculation. Since the N = 3 zig- 
zag tube essentially consists of three real-space chains, 
a spin gap may seem surprising. Indeed, it is knownE2l 
that conventional Heisenberg spin ladders with an odd 
number of legs have gapless spin excitations, in contrast 
to the spin-gapped even leg ladders. An important dis- 
tinction, however, is the unusual topology of the strong 
coupling Heisenberg model for the zig-zag tube, as de- 
picted in Fig. 4b. In the anisotropic limit, J± <C J, spins 
on such a "herringbone" lattice will indeed form local 
singlets across the vertical bonds, with a spin-gap. 

Returning to the armchair tubes, since the effective 
interaction strength un = u/N, one expects correlation 
effects to be weak for large N. Indeed, as noted above, 
the gaps in the undoped case become exponentially small 
for un < t, and the scale for superconductivity will like- 
wise be small, indicating the desirability of reducing the 
nanotube size in experiments. 

Even for larger N (the weak-coupling limit), how- 
ever, an interesting observable consequence of interaction 
physics should remain in the high-temperature resistiv- 
ity. Indeed, Umklapp scattering leads to an intrinsic con- 
tribution to the scattering rate which in weak coupling 
varies linearly in T for T > A c in ld!El This is a dramatic 
enhancement over the conventional Fermi liquid T 2 resis- 
tivity, and can also be much larger than scattering due 
to (3d) phonons, which vanishes at least as fast as T 3 . 
We now proceed to obtain a quantitative estimate of this 
effect. 

The effective low-energy Hamiltonian for the (N, N) 
nanotube (and for the two-chain Hubbard model) con- 
sists of right and left moving electrons in the two bands: 

H = ^2 / dx[i()\ ta ivd x 'ip Ra - ipl a ivd x ip La ], (7) 

a=l,2 J 

where we have suppressed the spin label. Since the equiv- 
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alent 2-chain Hubbard model is at half-filling, the pres- 
ence of the Hubbard interaction, un, introduces Umk- 
lapp scattering, as well as numerous momentum conserv- 
ing four-fermion interactions. The three Umklapp inter- 
actions, which scatter two right moving electrons into two 
left movers, take the form 



H v = u N a Q / dx^L^Ii ipmipin + h.c] 



(8) 



where we have now suppressed the band index. The scat- 
tering rate from Umklapp scattering can be extracted 
from the imaginary part of the electron's self energy, 
T(u), T) = ImT,(uj, T, k F ). To lowest order there is a sin- 
gle diagram for each of the three Umklapp interactions, 
which give identical contributions. One finds: 



F{u,T) = —(u N a /v) 2 TT(Lo/2t), 
an 



(9) 



where r(A") is a scaling function which approaches 1 as 
X — > 0, and varies as \X\ for large X. If we ignore vertex 
corrections, the Kubo formula for the Id conductivity 
can be expressed in terms of T as, 
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duj j-duf) 
2tt T(uj,T)' 



(10) 



where / = (e' 3 " + 1) 1 is the Fermi function. The result- 
ing Id resistivity is 



c h 

p(T) = — ^(u N a /hv) 2 (T/hv), 

167T e z 

with c a dimensionless constant of order one. 
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(11) 
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with metallic behavior. At lower temperatures the re- 
sistivity appears to saturate, perhaps turning up slightly 
but showing no compelling sign of a sizeable charge gap. 

A comparison with these results can be made by con- 
verting Eq. ( [TT| ) to the 3d resistivity pzd ~ pD 2 , where 
D is the nanotube diameter. This gives the rough es- 
timate p3d ~ 2(u/t) 2 (T/t)p£lcm. Notice that the nan- 
otube size N has dropped out. In the experiments, 
dpexp/dT w 10~ 2 pflcm / K . To account for this mag- 
nitude one would need a rather large bare Hubbard in- 
teraction, u/t ~ 10, perhaps not unreasonable given the 
neglect of long-ranged Coulomb forces in our simple Hub- 
bard treatment. The finite residual resistivity as T — > is 
presumably due to disorder. For example, local kinks or 
other defects in the rope packing would naturally lead to 
a temperature-independent additive contribution to the 
resistivity (see Fig. 5). However, other effects such as 3d 
crossover may also play a role at low temperatures. 
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FIG. 5. Resistivity of an ideal armchair tube (schematic). 

One of the most promising recent experiments! stud- 
ied single wall carbon nanotubes packed together into 
a triangular lattice to form crystalline "ropes". These 
ropes have diameters of 50 to 200 Angstroms, are tens to 
hundreds of microns long, and are believed to be predom- 
inantly composed of (10, 10) armchair tubes. Transport 
data on a single rope reveals a resistivity increasing lin- 
early with temperature between 50 and 300-ftT, consistent 
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